VISCOUS REGULARIZATION OF THE EULER EQUATIONS AND 

ENTROPY PRINCIPLES* 

JEAN-LUC GUERMONDt AND BOJAN POPOV* 

Abstract. This paper investigates a general class of viscous regularizations of the compressible 
Euler equations. A unique regularization is identified that is compatible with all the generalized 
- - - entropies a la Harten ct al. llOll and satisfies the minimum entropy principle. A connection with a 

^ recently proposed phenomenological model by Brenner [J] is made. 
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1 . Introduction. Proving positivity of the density and internal energy and prov- 
I ing a minimum principle on the specific entropy of numerical approximations of the 

compressible Euler equations is a challenging task that has so far been achieved for 
very few numerical schemes on arbitrary meshes in two and higher space dimensions. 
The Godunov scheme (Godunov 0) and some variants of the LaJ^ scheme (Lax 13 1) 



are known to satisfy all these properties, (see Einfcldt ct al. [2] for the Godunov 
scheme, Perthame and Shu [isl Appendix] for the explicit Lax algorithm, and Tang 
I and Xu [l^ for the implicit version of the Lax algorithm). The argumentation for the 

Godunov scheme relies on the fact that Riemann problems are solved exactly at each 
time step and averaging Riemann solutions preserves the above mentioned properties. 
None of the above arguments can be readily extended to central high-order schemes 
^ I and more generally to schemes that are based on Galcrkin approximations. One way 

\Q . to address this issue consists of using the standard parabolic regularization of the Eu- 

' ler equations to construct a scheme for which the vanishing viscosity is proportional 

jn , to the mesh size. The problem with this approach is that the regularization acts on 

the conserved variables which are the density, momentum, and total energy. Since 
. the momentum and total energy arc not Galilean invariant, a change of reference 

frame by translation and/or rotation changes the regularization. A way out of this 
dilemma consists of considering the Navier-Stokes regularization as a starting point to 
construct a numerical method. However, one then encounters two serious difhculties. 
The first one is that the Navier-Stokes equations do not include any regularization 
in the continuity equation, which is inconsistent with most numerical discretizations. 
^ . The second one is that whereas it is known that the Euler equations satisfy a mini- 

mum entropy principle on the specific entropy (sec e.g., Tadmor [l8|). it is also known 
that the Navier-Stokes equations violate this minimum principle if the thermal dif- 
fusivity is nonzero, see e.g., Serre [lil Thm 8.2.3]. These two observations make the 
Navier-Stokes regularization inconvenient for numerical purposes. One is then lead to 
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ponder on the following question: Is it possible to find a regularization of the Euler 
equations that is Galilean invariant, ensures positivity of the density and internal 
energy, satisfies a minimum entropy principle, and is compatible with a large class of 
entropies inequalities? The objective of this paper is to answer to this question. 

The paper is organized as follows. The parabolic and the Navier-Stokes regu- 
larizations and their apparent shortcomings mentioned above are discussed in fj2l A 
general family of regularizations is introduced and investigated in f|3] and 21 The 
minimum entropy principle is investigated in and the compatibility with entropy 
inequalities is studied in 21 The key result of this paper is Theorem 14.11 only one 
regularization technique satisfies the minimum entropy principle and is compatible 
with all the generalized entropies of Harten et al. . This formulation is compared 
in SjSlwith a reformulation of the Navicr-Stokcs equations proposed by Brenner [l| that 
is based on heuristic arguments. A striking observation is that by distinguishing the 
so-called mass and volume velocities, it is possible to re- write the proposed regulariza- 
tion into a form similar to that of the Navier- Stokes equations with rotation invariant 
viscous fiuxes. This way of looking at the regularization reconciles the parabolic and 
Navier-Stokes regularizations and shows that they are two faces of the same coin. The 
key results of the paper are summarized in ^5.31 Standard identities and inequalities 
from thermodynamics that are used in this paper are collected in Appendix [XI 

2. Standard regularizations. We review in this section some well-known reg- 
ularization techniques and discuss the pros and cons thereof. 

2.1. Statement of the problem. Consider the compressible Euler equations 
in conservative form in R^, 



where the dependent variables are the density, p, the momentum, m and the total 
energy, E. We adopt the usual convention that for any vectors a, 6, with entries 
{ai}i=i,...,d, {bi}i=i,...,d, the following holds: (a (g) b)ij = a^bj and V-a = dx^aj, 
{'^a)ij = dxiCLj- Moreover, for any order 2 tensors g, Ih, with entries 
{h^j}^.j=l....^d, we define (V-g)j = dx.gij, a-V = aidx,, (g-a)-t = .g-tjaj, g:[h = gijhij 
where repeated indices are summed from 1 to d. 

The pressure, p, is given by the equation of state which we assume to derive from 
a specific entropy, s(p, e), through the thermodynamics identity: 



9tP + V-m = 0, 
dtm + V-(i6 ® m) + Vp = 0, 
dtE + ^-{u{E + p)) = 0, 
p(a;,0) = po, m(a;,0) = mo. 



E{x,Q) = Eo, 



(2.1) 
(2.2) 
(2.3) 
(2.4) 



T ds := de + pdr, 



(2.5) 



where t :^ p ^, e := p ^E — ^u^ is the specific internal energy, u :~ p is the 

velocity of the fluid particles. For instance it is common to take s = log(e"'-i p~^) 
for an polytropic ideal gas. Using the notation Sg := || and Sp := this definition 
implies that 



Se-=T-\ Sp:=-pT-^p-\ (2.6) 



The equation of state takes the form p : 




or 



pSe + p^Sp = 0. 



(2.7) 
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The key structural assumption is that — s is strictly convex with respect to t :— p^^ 
and e. Upon introducing (T(r, e) := s{p,e), the convexity hypothesis is equivalent 
to assuming that ctt-t ^ 0, a^e < 0, and a-rr'^ee ~ f^e — (s'^c e.g., Godlewski and 
Raviart [^). This in turn implies that 

dpip'sp) < 0, See < 0, < dpip^Sp)See - P^ S%, (2.8) 

or equivalently that the following matrix 

Y^.^ (P^^dpip^sp) pspe\ ^2 9) 

\^ pSpe ps^e / ' 

is negative definite. In the rest of the paper we assume that (|2.8p holds and the 
temperature be positive 

0<Se. (2.10) 

Remark 2.1. Note in passing that contrary to what is sometimes done in the 
literature, we do not assume that the pressure be positive, which requires Sp < (see 
e.g., Godlewski and Raviart j^, p. 99], Harten et al. (2.3)]). For instance, the 
assumptions (|2.8|) and (|2.10|) hold for stiffened gases, but the quantity Sp can change 
sign. It is shown in the Appendix (see Remark lA.ll) that the convexity assumption 
(|2.8p and the positivity of the temperature (|2.10p are sufficient to prove that the Euler 
system is hyperbolic. This fact was first established by Godunov [8| in one dimension. 
It was established again in Friedrichs and Lax Q and Harten et al. [IB ]. 

The objective of the present paper is to introduce a viscous regularization of 
(|2.ip - (|2.4|) that is compatible with thermodynamics and that can serve as a reasonable 
starting point for numerical approximation. 

2.2. Monolithic parabolic regularization. A common regularization of (|2.ip 
for theoretical and numerical purposes consists of the following monolithic parabolic 
regularization: 

atp + V-m = eAp, (2.11) 

9tTn + V-(m ® m) + Vp = eAm, (2.12) 

dtE + Viu{E + p))^e/^E, (2.13) 

p(a;,0) = po, m(a;,0) = mo, E{x,0) ^ Eq, (2.14) 

where e is a small parameter. We call this regularization monolithic since no distinc- 
tion is made between the conserved quantities, i.e., the operator eA is applied blindly 
to all the conserved quantities. 

It can be shown that the Lax-Friedrichs scheme and its parabolic analog intro- 
duced in Perthame and Shu [l^] are approximations of (|2.1ip . For instance, consider- 
ing a nonhnear conservation equation dtU + \i'-F{U) = 0, where U is the dependent 
vector- valued variable in R™, the scheme introduced in Lax [3, p. 163] in one space 
dimension consists of considering 

U7^' - liu'Ui + uu) \kf{u-^^) F{UU)) (2.15) 
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where h is the mesh size, r is the time step, and A := Th~^. Assuming the flux 
F to be uniformly Lipschitz, to simplify, and upon introducing the maximum wave 
speed P := ||-F''||L~(Rm;R™ xr™) and the CFL number cfl := (3Th~^, (j2.15p is the 
centered second-order approximation of the following parabolic regularization of the 
conservation equation dtU + V-F{U) — eAU — 0, with the artificial viscosity e := 
i/iA~^ = ■^■^Ph. In other words, the Lax-Friedrichs scheme is a centered second- 
order approximation of (|2.1ip - (|2.14p with the numerical viscosity e = + 
c||L«:(R<ixiR+)j where c is the speed of sound. That the CFL number appears at the 
denominator of the artificial viscosity makes this scheme over-dissipative. It is often 
more appropriate to consider the following alternative 



which is also a centered second-order approximation of the parabolic regularization 
dtU + W-F{U) — eAU = with the viscosity ^f3h, which is more traditionally as- 
sociated with up-winding. This algorithm is often abusively referred to as the Lax- 
Friedrichs scheme. Both the above numerical schemes have interesting positivity and 
entropy properties, see e.g., Lax 12 1, Tadmor 3, 3], Perthame and Shu [ist . 



Despite its appealing mathematical properties, the above regularization is often 
criticized by physicists since it seemingly violates the Galilean and rotation invari- 
ance. It also dissipates the density, the momentum and the total energy, which 
seemingly are again aberrations from the physical point of view. When looking at 
()2.1ip - ()2.14|) . it is indeed difficult to see how this set of equations can be reconciled 
with the Navier-Stokes equations which arc usually viewed by physicists to be the 
acceptable regularization of the Euler equations. 

2.3. Navier-Stokes regularization. As mentioned above, a common "physi- 
cal" way to regularize the Euler system ()2.ip - (|2.4p consists of considering this system 
as the limit of the Navier-Stokes equations 

9tP + V-m = 0, (2.16) 

dtm + V-{u®m) + Vp-V-g = 0, (2.17) 

dtE + V-{u{E + p)) - V-{h + g-it) = 0, (2.18) 

p{x,0) = pQ, m{x,0) = mo, E{x,0) = Eq. (2.19) 



where g and h are the viscous and thermal fluxes. The most elementary model 
compatible with Galilean invariance consists of assuming that 

g = 2/iV"tt + AV-mD, h^KS/T. (2.20) 

where Vit := ^(Vm-I-(Vm)^), I is the identity matrix in R'', and T is the temperature, 
T := SjT^ . The viscosity ^ and the thermal diffusivity k are required to be non-negative 
by the Clausius-Duhem inequality, although these two parameters may depend on the 
state (p, e). 

We claim that (|2.16p - (|2.20p is not appropriate for numerical purposes and we 
identify at least two obstructions. The first problem is that the minimum entropy 
principle cannot be satisfied for general initial data if the thermal dissipation is not 
zero. More precisely, assuming k 7^ 0, for any r £ R, there exist initial data so that the 
set {s > r} is not positively invariant. Let us recall a simple proof of this statement 
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borrowed from Scrre [igI Thra 8.2.3]. The specific entropy for the Navier-Stokes 
system satisfies 

dts + M-Vs = 4;; (g:V"M + V-(kVT)) . (2.21) 

Assume that Mq := 'tioPo^^ constant. Assume also that the equation of state of the 
fluid is such that pe ^ 0, then one can use T and s as independent state variables 
since p^det ^ i?(p e j ) = ^{spSee — SeSpe) = Pe 7^ (scc ()A.6P ). One can then choose 
So with global minimum at and Tq so that ATo{0) < and Vro(O) = 0. Without 
loss of generality, we assume that k > in a neighborhood of 0. Then dts{0,0) = 
kPq^{0)ATo{0) < 0, thereby proving that {s > r} is not positively invariant for the 
regularized system (PTTB)) - (P^ . 

Another argument often invoked against the presence of thermal dissipation is 
that it is incompatible with symmetrization of the Navier-Stokes system when using 
the generalized entropies of Hartcn for polytropic ideal gases. (The function pf{s) 
is said to be a generalized entropy if f'j~^ - f" > 0, /' > and / e C^(IR;IR), 
see Harten Q.) It is proved in Hughes et al. [lli that the only generalized entropy 
that symmetrizes the Navier-Stokes system p.l6p ~ (|2.20p is the trivial one ps when 
K ^ 0, see also Tadmor [l^ (2.11) and Remark 2, page 460]. Note though that 
symmetrization of the viscous fluxes is not necessary to prove entropy dissipation. 
It is nevertheless true that the Navier-Stokes system with k 7^ docs not admit a 
generalized entropy inequality if /"(s) 7^ 0, and this fact is a consequence of the 
following quadratic form not being non-negative: f'{s)X^ — f"{s)XY, {X,Y) G R^. 
Symmetry of the viscous flux is not a necessary condition for entropy dissipation, see 
e.g., Serre [ll §1.1]. 

The above two arguments seem to imply that one should take k = if one wants 
to use the Navier-Stokes system as a numerical device that regularizes the Euler 
equations, satisfies the minimum entropy principle, and satisfies entropy inequalities. 
In that case, one then faces a serious obstruction when solving for contact waves. For 
instance assuming that the initial data, po, mo, Eq are such that the exact velocity is 
constant in time and space, say u = Pe^, the problem ()2.16|) - ()2.19p reduces to solving 
two linear transport equations 

dtp + Pd^P = Q, p(-,0) = po, (2.22) 
dtE + pd^E^O, E(-,0)^Eo. (2.23) 

Note that u being constant implies that the pressure gradient is zero. The exact 
solution is p{x, t) = po{x — Ptex). To make things a little bit more interesting assume 
that pq is piecewise constant, say ^0(2;) = 1 if a; < and ^0(2;) = 2 if a; > 0. In the 
absence of some sort of regularization, the above two linear transport equations are 
difficult to solve numerically. Except for the method of characteristics and Lagrangian 
based techniques, we are not aware of any numerical methods that can solve these 
equations without resorting to some kind of viscous regularization. 

In conclusion, if positivity of the density, the minimum entropy principle and 
a reasonable approximation of contact discontinuities is desired, the Navier-Stokes 
regularization does not seem to be appropriate to regularize p.l[) - ()2.4p . whether k is 
zero or not. 

3. General regularization. We investigate in this section the properties of a 
class of rcgularizations that we expect to be as general as possible. More precisely. 
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let us consider the following general regularization for the Euler system: 

+ V-m- V-/ = 0, (3.1) 
9fm + V-(tt(g)m) + Vp- V-g = 0, (3.2) 
dtE + Viu{E +p)) - Vih + g-u) = 0, (3.3) 

where for the time being we let the fluxes /, g, and h to be as general as possible. 
A theory of viscous regularization for general nonlinear hyperbolic system has been 
developed in Serre and Serre [l^ Chap 6] . This theory identifies classes of entropy- 
dissipative viscous regularizations and establishes short term existence results. Our 
objective in this paper is more restrictive. We want to construct the fluxes /, g, and 
h so that p.ip - p.3p gives a positive density, gives a minimum principle on the specific 
entropy, and is compatible with a large class of entropies. (Note in passing that the 
positivity of the internal energy will be a consequence of the positivity of the density 
and the minimum entropy principle.) In the rest of the paper, we are going to work 
under the assumption that (|3.ip - p.3p has a smooth solution. 

3.1. Positivity of the density. Let us now choose the flux / so that it regular- 
izes the mass conservation equation. From the theory of second-order elliptic equation 
we conjecture that a(p, e)V/9 should be appropriate, where a{p, e) is a smooth positive 
function of p and e. In particular, it is reasonable to expect that the following choice 
implies positivity of the density: 

a{p,e) = x{p,e)(p'{p), (3.4) 

where x is a smooth positive function of p and e and </? is a strictly increasing function. 
This definition gives / = x(yO, e)V(p(p). This regularization is at least compatible with 
the positive density principle as stated in the following. 

Lemma 3.1 (Positive Density Principle). Let f = a{p,e)Wp in (|3.ip . with a G 
L°°([R2;[R) and inf(5_^)gR2 a(C, r?) > 0. Assume that u and V-tt G L°°(R''xlR+; R). 
Assume also that there are constant states at infinity p°° , u°° , so that the supports of 
p{-, •) — p°° and u{-, •) — are compact in R'^x{0,t), for any t > 0. Assume finally 
that po — Pea G L^(IR'^; R). Then the solution of (|3.ip is such that 

ess inf p(a;,i) > 0, Vt > 0. (3.5) 

Proof. Owing to the assumed regularity of u and po, the theory of parabolic equa- 
tions implies that there is a unique solution to (|3.ip such that p—poo & L°°{{0, oo); i^(IR''))n| 
L^{{0,oo);H\R'^)) and dtp e L^iiO,oo); H-^iR'^)), see e.g., Evans d p.356]. 

Let e > and let h^^x) be a smooth concave function that approximates min(x, 0) 
uniformly over R; say there is c > so that sup^gj^ |/ie(s)— min(s, 0)\ + \hf{s) — sh'^{s)\ < 
ce and h" < 0. Let i > be some fixed time. Let B{0, R) be the ball centered at 
of radius R such that the supports of p(-, r) — p°° and u{-, r) — u°° are in i?(0, R) 
for all r G [0,t]. Let ^ be a regularized characteristic function with the following 
properties: x\b(o.b.) = 1 and x|Rrf\_B(o,i?+i) = 0. Multiplying the weak form of p.ip 
by the legitimate test function xh[{p) we obtain 

/ ({dth,{p) + uVKXp) + pK{pW-u)x{x) + aVp^ixKip))) dx = 0. 
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Using that the properties of x, we simplify the above equation as foUows: 

/ (dMp) + uVh,{p) + pK{p)\/-u + h'^{p)a\yp\^) dx = 
^ (dtK{p) + V\K{p)u) + {ph!Xp) - h,ip))V-u + /i'/(p)a|Vpp) da: = 



Now, we integrate over time and, owing to the assumptions regarding the behavior of 
u, p and a, we obtain 



/ h,{p{x,t))dx>- f f \ph'^{p)~h,{p)\\V-u\dxdt 

jRd Jo JR^ 



h^{pa{x)) 



> -ce+ / KXpo{x)). 



We can now pass to the hmit on e using the Lebesgue dominated eonvergence and we 
obtain min(p(a3, t), 0) > 0. The resuh follows readily. □ 

3.2. Minimum entropy principle. We now investigate under which conditions 
on the fluxes /, g and h, a minimum principle on the speeifie entropy holds. In order 
to account for impact of the viscous part in the mass conservation, we change the 
notation of the various viscous fluxes as stated in the following lemma. 

Lemma 3.2. Setting g = G + / ® m, and h = I ^u^ f , the specific entropy for 
the system (|3.ip - (|3.3p satisfies 

p[dts + M-Vs) + V-((ese — psp)f — Sel) — f-\i'{ese — psp) + Z-Vsc — SeG: Vu = 0. (3.6) 
Proof. We re-write p.ip - p.3p in non-conservative form as follows: 



dtp + u-Vp + pV-u - V f = 0, 

p{dtu + u-Vu) + uV f + Vp- V-g = 0, 

p{dt£ + M-V£) + £V-/ + V-(mp) - V-(/i + Q-u) = 0. 

where we have defined £ ~ p^^E. Then we obtain the equation controlling the internal 
energy, e = £ — ^u^, by multiplying the momentum equation by u and subtracting 
the result from the total energy equation: 

p{dte + M-Ve) + (e - ^u'^)\7-f + pV -u - V h - g:VM = 0, 

The key to obtain the equation that controls the entropy is to multiply the mass 
conservation by pSp, multiply the internal energy balance by Se, and add the two 
resulting equations. This linear combination is motivated by the following observation 
daS = SpdaP + SedaC which holds for any independent variable a € {t,x}. We then 
obtain 

p{dts + M-Vs) + Se(e - 5it^)V-/ + {pse + p^Sp)V-U 

- Se{V-h + g-Vtt) - pspV-f = 

The definition of the pressure implies that the quantity pse + p^Sp is zero, see (|2.7p . 
This simplification yields 

p{dts + tt-Vs) + (ese — psp)y-f — Se(g:Vtt) — Se^U^'^-f — Se^h = 0. 
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Wc now regroup the terms 

p{dts + u-Ws) + (ese - pSp)V-f - SeV-(/i + \v?f) - .Se(g:VM -{f® u):Vu) ^ 0, 

and conclude by using the definitions g — G(V*m) + f ®u and h = I ^u^f- D 

From now on we assume that the following structure holds for the viscous fluxes 
introduced in (|XT1) - (P3|) : 

g^G{V'u) + f®u, h^l-^u^f, G(V"m):Vm>0. (3.7) 

We also assume that / has the following form: 

/ = a(p,e)Vp a(p,e)>0, (3.8) 

and I is defined so that 

I ^ s^^{ese — psp)f + d{p,e)ps'^^\/s, d{p,e)>Q. (3.9) 

Remark 3.1. The conditions G(V*it):VM > 0, a{p,e) > 0, and d{p,e) > are 
essential to establish the minimum principle on the specific entropy and the entropy 
inequalities (see Theorem 13.41 and Theorem 14. ip . 

Remark 3.2. The structural assumption I ~ s'^^{eSe — psp)f + c?(p, e)psj^ Vs is 
crucial. This condition is equivalent to assuming that the conservative term in (|3.6p 
is of the following form: (esg — psp)f — s^l = —V-{dpVs). The definition of I makes 
sense since thermodynamics requires that Sf. = T^^ > 0, (see (|2.10p ). Note that given 
p.Sp the following alternative forms hold I = {d — a)pSpS^^W p + aeV p + dpVe^ or 
l = {a- d){pp-^ + e)Vp + dV{pe). 

Let us define the quantity 

J := -/•V(e.Se - psp) + l-\/se + aVp-Vs (3.10) 

which is a quadratic form with respect to Vp and Ve and whose coefficients depend 
on p, e, a{p,e), c{p,e), and d{p,e). 

Let Bd be the dxd identity matrix. For any symmetric 2x2 block matrix N 

^^(niih n,,^A denote := f"" "^^V 

Given row vectors X,Y E R'*, the quadratic form {X ,Y)-\t4-{X ,Y)'^ , generated by 
the 2x2 block matrix IM, is negative semi-definite if and only if N2 is negative semi- 
definite, i.e., 7122 < and det(N2) < 0. 

Lemma 3.3. Assume that (j3.8p - (|3.9p hold. The quadratic form J is negative 
semi-definite if and only if 

addet(S) -jid- afp'^slpl > 0. (3.11) 

Moreover, let X £ R such that d{l + X) = a, then 

J + Xd^Vse-Vs <0. (3.12) 



The inequality (j3.12p becomes strict if a > and d> 0. 
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Proof. Using the definition of we re- write J in the following form: 

J = —aseVpVe — aeVpVse + asp\V + apV pV Sp + aeVpVse — apspV pVse 

+ dps^^S/seispS/ p + SeVe) + aS/ p[sp\/ p + SgVe). 

This expression can be further simplified as follows: 

J = 2aSp|Vpp + ap\/p{spp\/p + Spe\7e) 

+ (d- a)pSpS^^Vp{SpeVp + SeeVe) + dpS/e{speVp + SeeVe) 

= (Vp, Ve)^N(Vp, Ve), 
where the matrix N is defined by 

7^11 = (rf - a)psps^'^spe + ap^^dpip'^Sp), 

2ni2 = (rf - a)pSpS^^See + (rf + a)pspe, 

Let us define the 2x2 block matrix Q obtained by setting a = and = 1 in N: 

911 = PSpS'^Spe, qi2 = pSpS'l^^See + pSpe, ^22 = PSee- 

Notice that this definition implies that the quadratic form induced by Q is 
(Vp, Ve)-Q-(Vp, Ve)^ = — Vs^-Vs. 

Now let us consider the following 2x2 block matrix IM = N + AdQ where A S R. Let 
us set d' ^ d{l + X) and observe that 

"T-ii {d' - a)psps'^^spe + ap~'^dp{p'^Sp), 
2mi2 = {d' - a)psps'^^s^e + {d' + a)pspe, 

m22 = d'pSee- 

Observe finally that J + Ad^Vsg-Vs = (Vp, Ve)-trl-(Vp, Ve)^. 

To have a negative semi-definite form we need TO22 = d' ps^e < 0, which means 
< d' since See < owing to the convexity assumption (|2.8p . We also need det(Ill2) 
to be non-negative, 

det(M2) = {{d' - a)psps~^spe + ap^^ dp{p^ Sp))d' pSee 

- \ {{d' - a)pSpSe^See + [d' + a)pSpeY 
= ad' {dp{p^Sp)See - P^S^e) - ji^,' - af p'^ S~'^ {SeSpe - SpSeef- 

Now if we set A so that d' = d{l + X) = a, then det(IM2) is non-negative and d' = a > 0. 
Note in passing that upon setting A = 0, this computation shows that J < if and 
only if ((31T1) holds. □ 

Remark 3.3. Note that we could avoid invoking the convexity of the entropy in the 
above argument by taking a = and A = — 1. This would however defeat the purpose 
of our enterprise whose primary goal is to find a nonzero viscous regularization of 
the mass conservation equation that ensures positivity of the density and is entropy 
compatible. 



N 
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Remark 3.4- Note that J < when a = d. 

Theorem 3.4 (Minimum Entropy Principle). Assume that pq and eo are con- 
stant outside some compact set. Assume also that p.7|) - (|3.8p - (|3.9p hold. Assume that 
the solution to (|3.ip - (j3.3[) is smooth, then the minimum entropy principle holds, 

inf s{x,t)> inf so{x), Vi > 0. 

Proof. We re- write p.6|) as foUows: 

p{dts + it-Vs) + V-((eSe — psp)f — sj) — f-\/{eSe — psp) + Z-Vsg — SgGiVit = 0. 

Upon using ()3.9p we obtain 

p{dts + it-Vs) - V-(dpVs) - f-V{eSe - pSp) + l-Vs^ - SeG:Vu = 0. 

Let N := —f-V{eSe — pSp) + l-VSf.. Using definition p.lOp we have N = J — aVp-Vs 
and 

p{dts + it-Vs) - V-(dpVs) - aVp-Vs = - J + SeG:Vu > 0. (3.13) 

Owing to Lemma [3.31 there is A e IR so that J + Ad— Vsg-Vs < 0. Finally we have 
proved that 

p{dts + u-Vs) - V-(dpVs) - {aVp + \d—Vse)-Vs 

(3.14) 

= - J - A(i— Vse-Vs + SeGlVw > 0. 

Se 

By assumption all the fields are smooth and s is constant outside some compact 
set (i.e., p and e are constant outside some time-dependent compact set since the 
initial data are constant outside a compact set and the speed of propagation is finite) . 
For each time t, s reaches its minimum; let x^i^(t) be one point where the minimum 
of s is reached, then \/s{xnun{t),t) ~ and As(xmin(i), i) > 0. The equation (|3.14p 
implies that 

pdts{{xmin{t),t)) - dpAs{ 

which in turn implies that p9ts((xniin(^), ^)) > 0, and we conclude that the minimum 
entropy principle holds. □ 

Remark 3.5. Note that the condition p. lip is not required to hold for the mini- 
mum principle to hold. 

4. Entropy inequalities. We investigate in this section whether the rcgular- 
ization of the Euler equations p.ip - p.3p is compatible with some or all generalized 
entropy inequalities. 

4.1. Generalized entropies. Let us consider all the generalized entropy iden- 
tified in Harten et al. [l3l. A function pf{s) is called a generalized entropy if / is 
twice differentiable and 

/'(,s)>0, ns)c-' - ns) > 0, V(p,e)eR2+, (4.1) 

where Cp{p,e) = TdTs{p,T) is the specific heat at constant pressure. It is shown 
in Harten et al. [13] that — p/(s) is strictly convex if and only if (14.11) holds, i.e.. 
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(|4.ip characterizes the maximal set of admissible entropies for the compressible Euler 
equations that are of the form pf{s). 

Theorem 4.1 (Entropy Inequalities). Assume that ((X7)) - (P^ - (P^ hold. Any 
weak solution to the regularized system (|3.ip - (|3.3p satisfies the entropy inequality 

dtipfis)) + V-iupfis) - dpVfis) - a/(s)Vp) > 0. (4.2) 

for all generalized entropies pf{s) if and only if a ~ d. 
Proof Let us multiply ((3i:^ by /'(s), 

p{dtf{s) + M-V/(s)) - V-(dpV/(s)) + dp/"(s)|Vs|2 - a/'(s)Vp-Vs 

+ Jf'{s)^f{s)s,G:Wu. 

We now multiply the mass conservation equation p.ip by /(s) and wc add the result 
to the above equation: 

dtipfis)) + Vfupfis)) - V-(dpV/(s) + a/(s)Vp) 

+ dp/"(s)|V(s)|2 + Jf'is) = fis)s,G:Vu 

We now investigate the sign of the quantity c?p/"(s)|Vsp + Jf'is). 
Owing to (|4.ip . wc have 

dp/"(s)|Vs|2 + Jf'is) < (rfpc;i|Vs|2 + J)f'is). (4.3) 

We now need to determine the sign of the quadratic form in the right hand side of 
the above inequality: 

dpCp^lVsp + J = dpCp^\spVp + SeVef + J 

dpcp\sl\\/ p\^ + 2spS^Vp-Ve + s^lVep) + J = dp(Vp, Ve)-S-(Vp, Ve)^, 

where the coefficients of the 2x2 block matrix S are defined as follows: 

dsn = dcp^s^ + [id - a)spS^^Spe + ap^^dpip'^ Sp)) 

2dsi2 — 2dCp^SpSe + [id - a)SpS'^^See + id+ a)Spe) 
ds22 = diCp^sl + See), 

and can be re-written into the following form 

dsn = dicp^s^ + p^^dpip'^Sp)) + (rf - a)s^^ {spSpe - s^p^'^dpip^ Sp)) 

2dsi2 = 2diCp^SpSe + Spe) + (rf - a)s~^ (SpSee - SeSpe) 
ds22 = diCp^sl + See). 

Then upon setting a; = 1 — ^ wc infer that 

Sll = hii + Xp^^SePp, 2si2 = 2/li2 + Xp^^SePe, S22 = ^22 (4.4) 

where the 2x2 matrix IH2 is defined by 

t\2= ( ^^"^^^ ^ P^^dpip'^Sp) SpSeCp^ + SpA 
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is shown to be negative in Lemnia [A.3l In particular we have S22 = /i22 = sicp^+s^^ < 
owing to the inequality CpT^ > 1 established in (jA.12|) . As a result, the matrix S is 
negative semi-definite if and only if the determinant of S2 is non-negative, 

det(S2) = /lll/l22 + xh22P~^SePp - {hl2 + \xp~'^S^Pef 

= det(IH2) + xp^^Se{h22Pp - hi2Pe) - jx'^p'^slpl. 
According to Lemma lA.31 we have det(IH2) and h22Pp — hi2Pe ~ 0. This proves that 

det(S2) = -Ix'p-^slpl 

In conclusion, S is negative semi-definite if and only if a; = 0, ie a = d. 

The above argument shows that dp/"(s)|Vsp + J f'{s) < if a = d . This proves 
that all the generalized entropy inequalities are satisfied if a = d. 

If a 7^ d we consider generahzed entropies such that /"(s) = (1 — e)f'{s)cp{s, p), 
e G (0, 1) (it is always possible to solve this ODE for any fixed value of p). For this 
subclass of generalized entropies, we have 

dp/"(s)|Vsp + Jf'is) = ((1 - e)dpcp i|Vs|2 + J)f'{s). (4.5) 

From the proof of Theorem (|4.ip . wc know that the quadratic form dpc~^|Vsp + J = 
dp(Vp, Ve)-S(p, e)-(V/9, Ve)"^ is negative semi-definite if and and only a = d. Let 
(p*,e*) G be a pair of positive numbers so that a{p*,e*) ^ d{p*,e*). Since the 
quadratic form generated by S(p*, e*) is not negative semi-definite, there exists a pair 
of row vectors X,Y eR'' so that (X, Y)-S{p* ,e*)-{X , Y)"^ > 0. It is always possible 
to choose e small enough so that 

(X, Y)-S{p*,e*)iX, Y) - ed*p*{c;)-'\s;X + s^Yl^fis*) > 0. 

Now we define an initial state so that in the neighborhood of the origin we have 
the following data: mo = 0, ^0(3^) = P* + x-X, eo(a;) = e* + x Y. Notice that 
with this choice Vixq = 0, Vpo = ^ and Veo = Y; therefore dpo/"(so)|V(so)P -I- 
J(po, eo)/'(so) — /'(so)se(po, eo)G:VMo > 0, which proves that the entropy inequality 
is violated at the origin close to the initial time. In conclusion a = d is a necessary 
condition for all the generalized entropy inequalities to be satisfied. □ 

Remark J^.l. Upon re-defining the velocity u = u -\- {d — a)V logp, the entropy 
inequality (|4.2p can be re- written into the following form 

dt{pf{s)) + V-{upf{s)) - VidpVpf{s)) > 0. (4.6) 

Remark 4-2. Theorem l4 . 1 1 proves that the family of regularization such that a = d 
is the most robust in the sense that it is the most dissipative. This result suggests 
that the choice a ~ d may be a very good candidate to construct a robust first-order 
numerical method for solving the compressible Euler equations. 

Corollary 4.2. Let a be a real number, a < 1, and assume that p.7p - (|3.8p - 
(|3.9p hold. Any weak solution to the regularized system p.ip - p.3p satisfies the entropy 
inequality (|4.2p for all the generalized entropies pf{s) such that /' > and ac~^ f > 

f" i/2r-2A3 < 1-| < 2T+2A^ where T = {l-a)det(E)p^s-^p-^ and A = T{l+T). 
Proof. Wc proceed as in the proof of Theorem l4.1l where we replace IH by H" where 
is substituted by ac~^. Upon replacing by ac~^ in the proof of Lemma |A.3[ 
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wc infer that dct(IH^) = (1 - a)p-'^dct{j:) and s^ih^^Pp - h^-^p^) = (1 - a) p''^ dci{Y:,) . 
Then by defining S" as in (|4.4I) . where IH is substituted by H". we obtain 

det(S^) = (1 - a)p-^det{T.) + - a)det(I]) - \x^p^^slpl 

= p-\{l - a)det(S)(l +x)- \x^p-^slpl), 

where we defined a; = 1 — |. Then upon setting F = (1 — a)det{'ll) p^ s'^'^p'^'^ and 
A = r(l + r), we conclude that the matrix S" is negative definite if 

2r - 2A^ < 1 - ^ < 2r + 2A^ 
a 

which ends the proof. □ 

Corollary 4.3. Any weak solution to the regularized system p.ip - p.3p satisfies 
the entropy inequality (|4.2p for the physical entropy ps (i.e., /(s) = s) i/2r — 2A2 < 
1 - ^ < 2r + 2A3 where T = det{Y.)p^s-^p-^ and A = r(l + T). 

Proof Take a = in Corollary or use ([3JT|) . □ 

4.2. Ideal gas. Let us illustrate the above theory in the case of ideal gases, 
i.e., s — \og{e~ p~^) with 7 > 1. We have = 7(7 — l)e, Cp = 7(7 — 1)^^, 
det(E) = (7 - ly^e-^, f = aVp, and I = jde(^ - 1 + i)Vp + dpVe. The range for 
the ratio ad^^ for Corollary 14.31 to hold is 

^(1 - V7) < 1 - ^ < + V7). (4.7) 

7— 1 (17—1 

In particular the choice 1 — ^ = ^ is clearly in the admissible range for the physical 

entropy inequality. This particular choice is such that I = dpVe and / = d'^^^^V p , 
i.e., I does involve any mass dissipation. 

5. Discussion. We show in this section that the regularization proposed above 
is a bridge between the Navier-Stokes and parabolic regularizations of the Euler equa- 
tions that reconciles the two point of views. 

5.1. Parabolic regularization. The first natural question that comes to mind 
is how different is the general regularization p.ip - (j3.3p from other known regulariza- 
tions. In particular how does it differ from the parabolic regularization (|2.1ip - (|2.14p ? 
The answer is given by the following, somewhat a priori frustrating result: 

Proposition 5.1. The parabolic regularization (|2.11|) - (j2.13p is identical to (|3.ip - 
([3:3)) with dSJll-dSH) where a = d = e, G = epVu. 

Proof. The equality a ~ e comes from the identification / = eVp in the mass 
conservation equation in ()2.1ip and (jSTj. The identity eVm = eVp (E) u + epVu 
implies that upon setting g = f ®u + G with G = epVu^ the momentum conservation 
equations in (|2.12p and p.2p are identical. Upon observing that 

g-M = v?f + G-it = eu^Vp + ]^epVv? e^^pu"^ + ^"^/> 
we obtain that 

eS/E = eV(pe) + \7^epu'^ = e\/{pe) - + S'"- 

The energy equations in (|2.13p and p.3p are identical if one sets h = I — \u^f, with 
I ~ eV(pe), meaning d = e. U 
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Remark 5.1. Even when a ~ d, one important interest of the class of regulariza- 
tion (|3.1l) - (l3.3p . when compared to the monohthic parabolic regularization, is that 
it decouples the regularization on the velocity from that on the density and internal 
energy. In particular the regularization on the velocity can be made rotation invariant 
by making the tensor G a function of the symmetric gradient W^u. This decoupling 
was not a priori evident (at least to us) when looking at the monolithic parabolic 
regularization (f2lT|) - (f2TT3l) . 

5.2. Connection v^rith phenomenological models. When introducing the 
structural assumptions (|3.7p - p.9p in the balance equations (|3.ip - p.3p we obtain the 
following system: 

dtp + V-m-V-f = 0, (5.1) 
dtm + V-(m «) m) + Vp - V-(G(V"tt) + / ® it) = 0, (5.2) 
dtE + V-{u{E +p)) - V-{1 + ^u^f + GiV-'u)-u) =0, (5.3) 

When looking at (|5.ip - (|5.3p it is not immediately clear how this system can be rec- 
onciled either with the Navier-Stokes regularization or with any phenomenological 
modeling of dissipation. 

It is remarkable that this exercise can actually been done by introducing the 
quantity ~ u — p^^ f . The conservation equations then becomes 

atP + V-(M„p) = 0, (5.4) 
dtm + V-{um ®m) + Vp- V-(G(V'*tt)) = 0, (5.5) 
dtE + V-{u^E) ~ V-(i - ef) + V-((pD - G(V*m))-m) = 0, (5.6) 

with again m = pu and E = pe+ ^pu^. It is surprising that this system involves two 
velocities. It is also somewhat surprising to observe that the above system resembles 
the Navier-Stokes regularization. In particular if one sets a = d, the term I ~ ef 
becomes dpVe, which upon assuming de = CydT, reduces to d{p,e)pCy\7T, i.e., one 
obtains Fourier's law: I — ef = d(p, e)pCyVT. 

During the preparation of this paper, it has been brought to our attention that the 
regularization model that we propose above somewhat resembles, at least formally, 
a model of fluid dynamics of Brenner [ij (see e.g., equations (1) to (5) in Brenner 
[l[). The author has derived the above system of conservation equations (up to some 
non-essential disagreement on the term I — ef) by invoking theoretical arguments from 
Ottinger [l3| and phenomenological considerations. The mathematical properties of 
this system have been investigated thoroughly by Feireisl and Vasseur [J]. Brenner 
has been defending the idea that it makes phenomenological sense to distinguish the 
so-called mass velocity, u^, from the so-called volume velocity, u, since 2004 (or so). 
We do not want to enter this debate, but this idea seems to be supported by our 
mathematical derivation (|5.4p - (|5.6|) that did not invoke any had oc phenomenological 
assumption. Recall that our primal motivation in this project is to find a regulariza- 
tion of the compressible Eulcr equations that can serve as a good numerical device, 
and by being good we mean that the model must give positive density, positive in- 
ternal energy, a minimum entropy principle and be compatible with a large class of 
entropy inequalities. 
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5.3. Conclusions. Let us finally rephrase our findings. In its most general form, 
the regularized system (|5.4p - (|5.6p can be re-written as follows: 

5tp + V-(M™p) = 0, (5.7) 

dtm + V-(it„ ® m) + Vp - V-(G'(V'it)) = 0, (5.8) 

at£; + V-(it„£:)- V-q + V-((pD-G(V*tt))-it) =0 (5.9) 

M„i = M — a(p, e)V logp (5.10) 

q = (a - d)pVlogp + dpVe, a(p, e) > 0, d(p, e) > 0. (5.11) 

It is established in Lemma 13.11 that the definition of / = a(p, e)VjO is compatible 
with the positive density principle. The particular form of q in (|5.1ip results from 
the definition of Z, see l|3.9p . which is required for the minimum entropy principle to 
hold, as established in Theorem l3.4l It is finally proved in Theorem 14.1 1 that the most 
robust rogularization, i.e., that which is compatible with all the generalized entropy 
a la Harten et al. [lOj, corresponds to the choice a ^ d. Various relaxations of the 
constraint a = d are described in Corollary 14.21 and in Corollary 14.31 As observed 
in ij5.11 the parabolic regularization can be put into the form (|5.7p - (|5.1ip with the 
particular choice G = aVit, which is not rotation invariant and uses the same viscosity 
coefficient for all fields. 

Appendix A. Primer in thermodynamics. We collect in this appendix stan- 
dard results from thermodynamics that are used in the paper. 

A.l. Chain rule. Let $ : ^ («,/?) , — > $(a,/3) = (0(a, ^), ?/'(a, /?)) e be 
a C^-diffeomorphism. The following holds: 



(A.l) 



In particular if (j){a, (3) = a we have 



9a/3(«,V) = -|4r^' d^l3ia,^)= . ,1 (A.2) 

A.2. Speed of sound. The square of the speed of sound is defined to be 

^~dpp{p,s), (A.3) 

i.e., is the partial derivative of the pressure as a function of the density and the 
specific entropy. Using the chain rule, this definition is equivalent to 

dpp{p, s) = dpp{p, e) 4- dep{p, e)dpe{p, s), (A.4) 

and using (|A.2p with a = p, (3 = tjj = s, one obtains 

^Pp-—Pe{p,e). (A.5) 

Using the following representations of and p^: 

Pe= p'^S~'^{SpSee- SeSpe), P p = S~'^ [p^ S pS pe - S^d^p^ S p)) , (A. 6) 

the expression (jA.5|) also gives 

= p^S~^{2Se,SpSpe - S^p^^ d{p'^ S p) - S^Sge). (A. 7) 
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A. 3. Convexity of the entropy, det(E). Let us define the following matrix 

y 3 pe See J 

which, up to the p factor, is the Hessian of the entropy with respect to the variables 
(p^^,e). The convexity assumption on the entropy implies that See and p^^dp{p^Sp) 
are negative. We have the following characterization of the determinant of E. 

det(S) = slippTe - PeTp). (A.9) 

To prove this statement, we observe that the following holds owing to (|A.6P : 

^e^e ^ee^ ^e^P — ^pe: 

slPe = P^{SpSee - SeSpe) S^Pp = p^ {spSpe - SeP^"^ dp{p^ S p)) . 

The result is now evident. 

A. 4. Specific heat at constant pressure . The specific heat at constant pres- 
sure is defined to be Cp(p, e) = TdTs{T,p). 

Lemma A.l. The quantities det(S), and Cp are related by 

Cpdet(S) = slc^. (A.IO) 

Proof. Using the chain rule, we can re- write the above definition as follows: 

Cp{p,e) = s~'^{spPT{p,T) + SeeT{p,T)). 

The change of variable formula (lA.ll) with the convention {a — p, (5 — e) and (0 = 
p,ip = T) gives 

Pt{p,T) ^ , eT{p,T) ^ 

PpTe - PeTp PpTe - PeTp 

We then have the following expression for Cp 



-1 



{PpSe-PeSp) 
PpTe - PeTp ' 



(A.ll) 



Then using the expression of in (jA.5[) and the relation (jA.Qp . we arrive at the 
desired expression. □ 

Lemma A. 2. The following holds: 

CpT, > 1. (A.12) 

Proof. The definition of Cp implies that we need to estimate TsT{p,T)Te{p, e). 
The chain rule implies 

1 = Tse(p, e) = Tsp{p, T)pe{p, e) + Tst{p, T)Te{p, e). 

The result will be established if we can prove that Sp{p,T)pe{p,e) < 0. We now 
calculate Sp{p,T). The chain rule implies again that 

{sp{p,T)y^ =Ps{s,T) ^pp{p,e)ps{s,T) +pe{p,e)es{s,T). 
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Then using (jA.ip witli the convention [a — p, P = e) and {(j> ^ s , ip — T) gives 

Te — T 

Ps{s,T) = ^- — - — — , es(s,T)- 



sT— sT sT— sT 



This in turn implies that 

PpTe-peTp _ s^^det(i;) 



(sp(p,T))-i=p,,(,s,T) 



SpTe - SeTp p ^p, 



since SpT^ — SeTp = s^^(— SpSge + SeSpe) = —p '^Pe, whcrc wc uscd (|A.6p . In conclusion 
Sp{p,T)pe{p,e) = -sIpIp^^ det{j:y^ < 0, owing to and (|2.10p . which concludes 
the proof. □ 

Remark A.l. Note in passing that the convexity assumption p.Sp implies that 
Te > 0, which owing to (jA.12[) implies that Cp > 0. This in turn implies that > 
owing to (jA.lOp . i.e., the Euler system p.ip - (|2.4p is hyperbolic under the convexity 
assumption ()2.8p and the positivity assumption on the temperature ()2.10p . Positivity 
of the pressure is not needed to establish this fact. 

A. 5. Matrix IH2. Investigations on entropy inequalities involve the quadratic 
form induced by the matrix IH2 



lU _/^oCp^+p '^dp{p'^Sp) SpSpCp^+Spe 
ln9 — I _i 2 —1 

^p^e^p ^pe ^e^P ^ee 

Some key properties of this matrix are collected in the following lemma. 
Lemma A. 3. The following hold: 

(i) det(lH2) = 0. 

(ii) IH2 is negative semi- definite. 

(Hi) h22Pp - hi2Pe = 0- 

Proof ([H Using the expressions (|A.7p and (|A.9p for the speed of sound, c^, and 
det(S]), and the relation (jA.lOp . the determinant of IH2 is re- written as follows: 

det(H2) = (s^Cp^ + p-^dp{p^Sp)){slcp^ + See) - {SpSeCp^ + Sp^f 
= /?"^det(S) + Cp^(slp''^dp{p'^Sp) + s'^^S^e - 2SpSeSpe) 



p"^det(S) - Cp'cy 



This is essentially the result established in Harten et al. jlOl . p. 2126]. 

(jn]) Owing to the inequality 1 < CpTe established in (|A.12p . we infer that /122 = 
SeCp^ + See < 0, which together with proves statement dn]). 
Let us compute s~^{h22Pp — hi2Pe), 

Se'^{h22Pp - hi2Pe) = (Cp ^ - Te)pp - {s pS'^'^ c:^'^ - Tp)pe 

= PeTp -ppTe + Cp^S'^^{SePp - S pPe) ■ 

This proves that s~^{h22Pp — /ii2Pe) — owing to (jA.lip . □ 
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